Scaling behavior of domain walls at the T = ferromagnet to spin-glass transition 
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We study domain-wall excitations in two-dimensional random-bond Ising spin systems on a square 
lattice with side length L, subject to two different continuous disorder distributions. In both cases an 
adjustable parameter allows to tune the disorder so as to yield a transition from a spin-glass ordered 
ground state to a ferromagnetic groundstate. We formulate an auxiliary graph-theoretical problem 
in which domain walls are given by undirected shortest paths with possibly negative distances. Due 
to the details of the mapping, standard shortest-path algorithms cannot be applied. To solve such 
shortest-path problems we have to apply minimum-weight perfect-matching algorithms. We first 
locate the critical values of the disorder parameters, where the ferromagnet to spin-glass transition 
occurs for the two types of the disorder. For certain values of the disorder parameters close to the 
respective critical point, we investigate the system size dependence of the width of the the average 
domain- wall energy (~ L^) and the average domain- wall length (~ L'*' ). Performing a finite-size 
scaling analysis for systems with a side length up to L = 512, we find that both exponents remain 
constant in the spin-glass phase, i.e. 9 k, —0.28 and df ~ 1.275. This is consistent with conformal 
field theory, where it seems to be possible to relate the exponents from the analysis of Stochastic 
Loewner evolutions (SLEs) via di — 1 = 3/ [4(3 -f 5)]. Finally, we characterize the transition in terms of 
ferromagnetic clusters of spins that form, as one proceeds from spin-glass ordered to ferromagnetic 
ground states. 

PACS numbers: 75.50.Lk, 02.60.Pn, 75.40.Mg, 75.10.Nr 
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I. INTRODUCTION 

Ising spin glasses (ISGs) are among the most-basic 
models of disordered systems that allow for the study 
of phase transitions in the presence of quenched disor- 
der. ISGs are elaborately studied in statistical physics 
[U, 0, i, il and despite several decades of active re- 
search they attract a constant interest, challenging with 
still not well understood traits and unresolved questions. 
In the scope of this paper we investigate ground state 
(GS) spin configurations and minimum-energy domain- 
wall (MEDW) excitations in a 2d random-bond ISG. In 
brief, MEDWs are topological excitations that are in- 
duced by a change of the boundary conditions (BCs) from 
periodic to antiperiodic along one boundary of the sys- 
tem. In particular wc are interested in the scaling prop- 
erties of MEDWs close to the critical point at which the 
T = spin glass (SG) to ferromagnet (FM) transition oc- 
curs. From a phcnomcnological point of view, the physics 
of the SG ordered phase of ISGs with short ranged in- 
teractions, like the 2d model considered here, can be de- 
scribed in terms of the droplet scaling picture 0, H, 0] ■ 
Therein, the low-temperature behavior is dominated by 
droplet excitations, i.e. clusters of spins that are flipped 
relative to the GS spin configuration. Within the droplet 
picture, excitations like MEDWs posses an excitation en- 
ergy Ai? that scales with system size L as Ai? ~ L^, 
where 6 is referred to as stiffness exponent. The value of 6 
is assumed to be universal for all types of excitations and 
constant within the whole SG phase. Moreover, in a sta- 
tistical sense, they are self-similar fractals characterized 
by a fractal dimension c?f that is defined by the scaling 
of the average MEDW length as {I) ~ L"^^ . The advan- 



tage of working at zero temperature is that the GS prob- 
lem for the particular setup studied here can be solved 
by means of exact combinatorial-optimization algorithms 
d, [^, [13, [m, whose running time increases only poly- 
nomially with the system size. Hence, very large systems 
can be treated exactly, giving very precise and reliable 
estimates for the observables. For 2c? lattices, where the 
interaction strengths (bonds) between adjacent spins are 
drawn from a Gaussian distribution with zero mean and 
unit width, domain wall (DW) calculations using such al- 
gorithms resulted in the estimates 6* = -0.287(4) [T3l.[T3| 
and (if = 1.274(2) [31 . The negative value of the stiffness 
exponent indicates that the excitation energy required 
to introduce a MEDW gets negligibly small as L oo 
and thus, thermal fluctuations prevent a spin-glass or- 
dering for any non-zero temperature. The above value of 
the stiffness exponent was later on confirmed for contin- 
uous disorder distributions different from the Gaussian 
bond distribution [iGj . for droplet excitations respecting 
a Gaussian distribution of the bonds [l3, [l^, [3| and 
quite recently also for droplets within the ± J model j20| . 
Furthermore, recent studies suggested that MEDWs re- 
specting a Gaussian distribution of the bonds can be de- 
scribed by stochastic Loewner evolutions (SLEs) [2ll.[2^. 
SLEs are generated by a stochastic differential equation 
driven by a brownian motion. They describe the con- 
tinuum limit for various 2d random curves and their geo- 
metric properties relate to the statistics of several critical 
interfaces [23]. Within conformal field theory it further 
seems to be possible to relate the DW fractal dimen- 
sion to the stiffness exponent by means of the relation 
df — l = 3/[4(3-f6')], subsequently referred to as SLE scal- 
ing relation. For the pure spin glass, this is in agreement 



with the numerical estimates of 6 and df stated above. 

Here, we consider a random-bond Ising model that al- 
lows us to investigate the SG to FM transition at zero 
temperature, by tuning the mean value of the underly- 
ing disorder distribution. In a previous work, the related 
±J model was studied in 2d There exact matching 
algorithms to find ground states (GSs) where applied. It 
was found that, in the limit of large system sizes, the SG 
to FM transition occurs at a fraction pc = 0.103(1) of an- 
tifcrromagnctic bonds (—J) among ferromagnetic bonds 
(+J). Further, the critical exponents v and /3 that de- 
scribe the divergence of the correlation length and the 
order parameter, where found to be = 1.55(1) and 
/? = 0.09(1). Due to the discreteness of the distribution, 
the DWs are not unique and cannot be sampled in equi- 
librium for large systems. Hence the fractal dimension 
has not been determined in a precise way so far. 

To clarify whether the SLE scaling relation above holds 
within the whole spin-glass phase, wc use two different 
continuous distributions of the disorder, which allow us 
to calculate the fractal dimension df with high precision. 
For this purpose, we perform GS calculations by means of 
exact combinatorial-optimization algorithms and study 
the scaling behavior of MEDWs close to the critical point 
where the SG to FM transition occurs. At first, we per- 
form a finite-size scaling analysis for systems of moderate 
sizes [L < 64) to locate the critical points at which the 
transitions takes place. Then we perform additional sim- 
ulations for large systems (L<512) close to and directly 
at the critical points, to get a grip on the scaling behav- 
ior of the MEDWs. Finally, we characterize the tran- 
sition using a finite-size scaling analysis for the largest 
and second-largest ferromagnetic clusters of spins within 
the GS spin configurations. These clusters form as one 
proceeds from spin-glass ordered to ferromagnetic ground 
states. To summarize our results: we find that the SLE 
scaling relation holds in the SG phase up to a point very 
close to the respective critical points, but not right at the 
critical points. Moreover, MEDWs in the SG phase scale 
like self-similar fractals, while MEDWs in the ferromag- 
netic phase display a self-afhne scaling behavior. 

The paper is organized as follows. In section |IT] we 
introduce the model and describe the algorithmic tech- 
niques wc have used in order to obtain MEDWs. In sec- 
tion mil we present the results of our numerical simula- 
tions. We conclude with a summary in section [TVl 



II. MODEL AND METHOD 

We performed GS calculations for two-dimensional 
random-bond Ising spin systems with nearest-neighbor 
interactions. The respective model consists oi N = Lx L 
spins a ~ (cti, . . . , cttv) with at = ±1, located on the 
sites of a regular square lattice. The energy of a given 
spin configuration is measured by the Edwards- Anderson 
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FIG. 1: Domain wall samples for systems of side- length L = 
64 and different values of the disorder parameter p (Model I) . 
The samples are taken in the SG phase (p = 0.60), right at 
the critical point (p = 0.66) and in the FM phase (p = 0.72). 
Besides the system size L, the DW length I and its roughness 
h are illustrated. 

Hamiltonian 

Hi<y) -^2-^'^ ^''^1' (1) 

where the sum runs over all pairs of adjacent spins with 
periodic BCs in the a:-direction and free BCs in the y- 
dircction. Therein, the bonds Jij are quenched random 
variables drawn from a given disorder distribution. Sub- 
sequently, we distinguish two types of the bond disorder: 

(1) Model I, where one realization of the disorder con- 
sists of a random fraction p of ferromagnetic bonds and a 
fraction (1 — p) of bonds that are drawn from a Gaussian 
distribution with zero mean and unit variance, i.e. 

Pi(J) = (1-p) exp(-jV2)/\/2^ + p6{J-l). (2) 

There exists a critical value pc of the disorder param- 
eter that separates a spin-glass phase (p < pc) from a 
ferromagnetic phase (p > pc). As limiting cases we can 
identify the pure Ising SG at p = and the ordinary 
Ising ferromagnet at p = 1. A similar type of disorder 
was used earlier for Monte Carlo simulations that where 
carried out to study the FM to SG transition in 3d and to 
numerically verify the absence of an equilibrium "mixed" 
ferromagnetic-SG phase for the respective model [2^ . 

(2) Model II, where the bond strengths are drawn from 
a Gaussian distribution with mean pj and width crj, i.e. 

Pii(J) = exp(-(J - pjf/{2aj))/iV2^aj). (3) 

As a function of the reduced variable r = aj / fij we ex- 
pect to find a ferromagnetic phase (spin-glass phase) for 
r < Tc (r > Tc). An earlier DW renormalization-group 
study of small systems [25| supported by transfer-matrix 
calculations reported, amongst other things, a zero tem- 
perature FM to SG transition at Tc = 0.961(10) with 
J/ = 1.42(8). Further, for the pure SG (pj = 0), an ex- 
trapolation of the DW free energy to zero temperature 
resulted in a stiffness exponent 9 = —0.281(5). The frac- 
tal properties of the DWs were not studied in this work. 

In the above two models, the bonds are allowed to take 
cither sign, where a value Jij > signifies a ferromagnetic 
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coupling that prefers a parallel alignment of the coupled 
spins, while a value Jy < indicates an antiferromag- 
netic coupling in favor of antiparallel aligned spins. The 
competing nature of these interactions gives rise to frus- 
tration. A plaquette, i.e. an elementary square on the 
lattice, is said to be frustrated if it is bordered by an odd 
number of antiferromagnetic bonds. In effect, frustration 
rules out a GS in which all the bonds arc satisfied. 

Here, our intention is to get a grip on the geometric 
properties of minimum-energy DWs. These are topologi- 
cal excitations that are defined, for each realization of the 
bond disorder, relative to two spin configurations: Up, a 
GS spin configuration with respect to periodic BCs fur- 
ther characterized by the configurational energy Ep and 
(Tap, a GS respecting antiperiodic BCs characterized by 
the energy Eap. Antiperiodic BCs are realized by invert- 
ing the sign of all the bonds along one column in the 
x-direction. Comparing the orientation of the spins in 
the two GSs, one can distinguish two regions on the lat- 
tice: one where the orientation of a spin is the same in 
both GSs and another, where the orientation of a spin 
differs regarding the two GSs. Within these regions, the 
bonds between adjacent spins are either satisfied or bro- 
ken in both GSs likewise. Bonds that connect spins that 
belong to different regions on the lattice are satisfied in 
exactly one of the two GSs. The MEDW is the interface 
in between the two regions and as such, it runs perpen- 
dicular to the latter bonds. It has the property that its 
excitation energy 5E = Eap — Ep is minimal among all 
possible DWs that span the system in the direction with 
the free BCs. The basic observables related to a DW 
are its over all length its roughness h and its excita- 
tion energy 5E. MEDWs for three different values of the 
disorder parameter p introduced above (see Model I) are 
illustrated in figure [1] 

We now give a brief description of the algorithm that 
we used to determine the MEDWs. A more extensive de- 
scription of the individual steps of the algorithm can be 
found in [l5j . For a given realization of the bond disor- 
der, we first determine a GS spin configuration consistent 
with periodic BCs in the x-direction. Besides the mag- 
netization niL = IX^i'^il/^^ ^'^'i ^^'^ energy, this tells 
which bonds are satisfied/broken in the GS for that par- 
ticular disorder sample. For the 2d ISG on planar lattice 
graphs, i.e. when there are periodic BCs in at most one 
direction, exact GS spin configurations can be found in 
polynomial time. This is possible through a mapping to 
an appropriate minimum-weight perfect-matching prob- 
lem [lol . ITl| . a combinatorial-optimization problem 
known from computer science. Here, we state only the 
general idea of this method. For this mapping, the spin 
system needs to be represented by its frustrated plaque- 
ttes and paths connecting those pairwise, i.e. matching 
them. In doing so, individual path segments are con- 
fined to run perpendicular across bonds on the spin lat- 
tice. Those bonds that are crossed by path segments 
are not satisfied in the corresponding spin configuration. 
The weight of the matching is just the sum of the ab- 



solute values of all bond strengths that relate to unsat- 
isfied bonds. Hence, finding a minimum-weight perfect 
matching on the graph of frustrated plaquettes then cor- 
responds to finding a spin configuration with a minimal 
configurational energy, hence a GS. The use of this ap- 
proach permits the treatment of large systems, easily up 
to L = 512, on single processor systems. This GS spin 
configuration can further be used to set up a weighted 
dual of the spin lattice, whose weighted edges comprise 
all possible DW segments. The weighted dual is con- 
structed as follows: set up a new graph G ~ [V,E,lo), 
whose sites i £ V relate to the elementary plaquettes 
on the spin lattice. Its necessary to introduce two ex- 
tra sites that account for the free BCs. Two sites are 
joined by an undirected edge e £ E, if the corresponding 
plaquettes have a bond in common. For the weight as- 
signment on the dual, consider a bond on the spin lattice 
having a coupling strength Jy . If the bond is satisfied 
(broken) regarding the GS, the corresponding dual edge 
e gets a weight uj{e) ~ — 2|Jy| {uj{e) = -|-2|Jy|). The 
weighted dual now comprises all possible DW segments, 
where the weight of an edge is equal to the amount of 
energy that it would contribute to a DW. Every possi- 
ble DW links both extra sites on the dual, where the 
energy of a DW is the sum of the weights along the ac- 
cording lattice path. So as to have minimum energy, it is 
beneficial for a DW to include (avoid) edges with a neg- 
ative (positive) edgeweight. Consequently, a MEDW is a 
minimum- weight path on the dual that joins both extra 
sites. The dual G is an undirected graph that allows for 
negative edge weights and so as to construct minimum- 
weight paths on G, it requires matching techniques [2^. 
Therefore we need to map the dual to an auxiliary graph 
Ga and find a minimum-weight perfect matching on Ga 
which we finally can relate to a minimum-weight path on 
G. For each realization of the disorder, this procedure 
yields an explicit representation of the minimum-energy 
DW that we can easily probe for its geometric proper- 
tics. A more detailed description of the algorithm can be 
found in p^ . In the following we will use the procedure 
outlined above to investigate MEDWs for the random- 
bond Ising models introduced earlier. 



III. RESULTS 

So as to characterize the scaling behavior of MEDWs 
for the two disorder distributions introduced above, we 
first of all need to find the critical values pc (Model I) 
and Tc (Model II) of the disorder parameters at which 
the r = SG to FM transition takes place. Reliable es- 
timates for the location of the critical points can already 
be obtained from comparatively small system sizes, here 
we use L = 24,32,48,64. In general, one has to find a 
proper balance of system size and sample numbers that 
affect finite-size effects and statistical error, respectively 
[27| . Subsequently we can probe the asymptotic scaling 
behavior of the MEDWs at fixed values of the disorder 
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FIG. 2: Results of the FSS analysis for the Binder parame- 
ter bL associated with the magnetization, for different system 
sizes L. The main plot shows the unsealed data near the criti- 
cal point, while the inset illustrates the data collapse obtained 
for (a) Model I: pc = 0.660(1) and u = lA9{7), (b) Model II: 
rc = 0.970(2) and i/ = 1.49(4). 



parameters close to the critical points for large system 
sizes up to L = 512. 



Finite-size scaling analysis to characterize the 
T = spin glass to ferromagnet transition 



First, we will discuss the results for the Model I disor- 
der and afterwards report the results for the Model II dis- 
order more briefly. As pointed out above, at large values 
of p, there exists an ordered ferromagnetic phase, while 
for small values of p a spin-glass phase exists. Therefore, 



a proper order parameter to characterize the respective 
SG-FM transition is the magnetization rriL = \ J2i L"^ 
for a system of size L. In the following, we perform a 
finite-size scaling analysis (FSS) in order to locate the 
critical point pc and also estimate the critical exponents 
that describe the scaling behavior of the magnetization 
at criticality. The Binder parameter [2^ associated with 
the magnetization reads 
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A? 



(4) 



and is expected to scale as bL{p) ^ fi[{p — pc)L^/^], where 
/i is a size-independent function and v signifies the crit- 
ical exponent that describes the divergence of the corre- 
lation length as the critical point is approached. Here, 
we simulated square systems of size L = 24,32,48,64 
at various values of the disorder parameter p. Observ- 
ables are averaged over up to 3x10^ (2x10*) samples for 
the smallest (largest) systems and we utilized the data 
collapse anticipated by the scaling assumption above to 
obtain pc = 0.660(1) and v = 1.49(7) with a quality 
S = 1.25 of the data collapse [2§|, see figure [2ja). The 
value of the critical exponent v agrees within errorbars 
with the value ly = 1.42(8) obtained using a transfer- 
matrix approach [2^ . Note that both, the numerical 
values of pc and further agree with those that charac- 
terize the negative- weight percolation of loops and paths 
on 2c? lattices [sO], highlighting the close connection of 
the two optimization problems. The order parameter of 
the transition is expected to scale conform with the as- 
sumption rriLip)^ L~'^^'^ f2[{p — Pc)L^^''], h being a size- 
independent function, where the magnetization exponent 
f3 was obtained after fixing v and pc to the values stated 
above. The most satisfactory data collapse (5'= 1.83) was 
obtained using /3 = 0. 097(6), see figure[3] In general, the 
above scaling relation holds best near the critical point 
and one can expect that there are corrections to scal- 
ing off criticality. As a remedy, we restricted the latter 
scaling analysis to the interval [—0.5, +0.2], enclosing the 
critical point on the rescaled abscissa. Note that the val- 
ues for the exponents found here agree with those found 
from GS calculation for the ± J-model [3l| within the er- 
rorbars. Further, the exponents appear to be consistent 
with those that describe the paramagnet to ferromagnet 
transition for the random bond Ising model, regarding 
finite temperatures T < T* below the temperature T* 
that characterizes the multicritical Nishimori point [33 |. 
In the respective study, the exponents v = 1.50(4) and 
[3 — 0.095(5) where measured by means of monte carlo 
simulations for the random bond ± J Ising model on lat- 
tices with L < 64 at fixed T, while varying the fraction 
of ferromagnetic bonds on the lattice. As an alternative 
order parameter, we also studied the average path length 
{I) of the MEDWs, where we expect a scaling of the form 



(0^L'^V3[(P-Pc)i^/1 



(5) 



Therein, df signifies the fractal dimension of the DWs 
at the critical point and /a is another size-independent 
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FIG. 3: Results of the FSS for the average magnetization 
mhip) for different system sizes L for Model I disorder. The 
main plot shows the unsealed data near the critical point, 
while the inset illustrates the data collapse obtained for the 
parameters pc = 0. 660(1), 1^ = 1.49(7) and /3 = 0. 097(6). 



function. From a finite-size scaling analysis restricted 
to the interval [—0.75, +0.5] on the rescaled abscissa, we 
obtained df = 1.222(4) with a quality 5* = 1.33, see figure 
U (Note that for a more clear presentation, the argument 
along the abscissa in figure [4] reads \p — pc\L^^^)- For 
the somewhat larger interval [— 1,+0.5] we found = 
1.223(4) with 5"= 1.40 in agreement with the above value. 
Since we expect the average MEDW length at p = to 
scale as (Z) ^L'*'' (here, p = corresponds to the pure spin 
glass studied in [131), where df = 1.274(2), we can further 
estimate the asymptotic behavior fz{x) ~ ^v(ds-df) 
the scaling function in Eq. ^ as x— >— oo. This can be 
seen from the top branch in figure 51 where the function 
fz{x)'^ x^ '^^^^^ is shown as solid line and agrees well with 
the data. Note that via Eq. ([5|) the DWs at pc exhibit 
the fractal dimension dj, while for all values p < pc, 
the fractal dimension is given by df. Hence, the scaling 
ansatz Eq. ([5]) is based on the assumption that behavior 
in the SG phase is universal, which is tested below for 
much larger systems explicitly. 

For the "ferromagnetic" branch {x —>■ +oo), a simi- 
lar consideration yields the asymptotic scaling f3{x) ~ 
2,-0.33(1)^ indicated as a dashed line in figure [H 

Further, we found that the probability Pl{p) that 
the MEDW roughness is equal to L scales as Pl{p) = 
/4[(p — Pc)L^^^]^ shown in the inset of figure SI In the 
ferromagnetic phase the value of Pl tends towards zero 
and in the spin-glass phase it saturates around Pl ~ 0.12. 
Hence, as pointed out in [2^ . an asymptotic nonzero 
probability that the MEDW roughness is 0{L) can be 
used as an order parameter to detect the SG phase. 

Regarding Model II, we simulated systems of size L — 
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FIG. 4: Results of the FSS analysis for Model I disorder. The 
main plot illustrates the FSS of the average MEDW length 
If) for different system sizes L, where the best data collapse is 
obtained for the parameters pc = 0. 660(1), 1^ = 1.49(7) anddf = 
1.222(4). The solid an dashed lines illustrate the asymptotic 
scaling behavior of both branches as described in the text. 
The inset shows the scaling of the probability Pl (p) that the 
roughness of the MEDW is equal to the system size L. 



24, 32, 48, 64 at different values of the disorder parameter 
r. Here, we fixed the width of the disorder distribution 
to the value aj = l and we vary only its mean p,j . Ob- 
servables are averaged over up to 3 x 10^ (2 x 10"*) sam- 
ples for the smallest (largest) systems and we utilized the 
data collapse anticipated by the scaling assumptions for 
the Binder parameter (see figure HUb)) and the magne- 
tization (not shown) to obtain the values Tc = 0.970(2), 
j/= 1.49(4) (5"= 1.0) and/? = 0.09(l) (S' = 0.46). Note that 
the numerical values of and v agree within errorbars 
with the values = 0.961(10) and v = 1.42(8) obtained 
using a transfer-matrix approach [2^ . The scaling of 
the average MEDW length here yields a numerical value 
df = 1.249(5) (5"= 1.99) which can only be considered as 
an effect of the finite system size, see the discussion be- 
low. Further, the probability that the MEDW roughness 
equals L tends towards Pl ~ 0.12 in the SG phase, in 
agreement with the above results. 



B. Scaling behavior at fixed values of p and r 

We have carried out further simulations at a couple of 
selected values of p and r, see tables U and [Hi in order 
to probe the asymptotic scaling behavior of MEDWs re- 
garding the two disorder distributions introduced above. 
We therefore considered systems of size up to L = 512 
with lO'^ realizations of the disorder. In particular, we 
are interested in the asymptotic scaling behavior of the 
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average MEDW length (/) with respect to the system size 
L, defining the DW fractal dimension d{ via {I) ^L'^^ . We 
further study the scaling of the average MEDW rough- 
ness (/i), i.e. the extension of the lattice path in the 
direction of the periodic BCs, that defines the rough- 
ness exponent dr by means of {h) ^ L'^' . Both these 
obscrvables relate only to the geometric properties of 
the MEDW, see figure [1] Finally, wc investigate the 
size scaling of the mean AE = {\SE\) ~ L^'^ and width 
a{SE) = ^/{SE^) - {SEy - L"^ of the distribution of 
MEDW excitation energies. 

Again, wc first discuss the results for the Model I dis- 
order and afterwards state the results for the Model II 
disorder more briefly. The asymptotic scaling behavior 
of the average DW length allows one to obtain the frac- 
tal dimension by using a direct fit to the power law data 
over the entire range of system sizes L. A reliable and 
more systematic alternative is to investigate a sequence 
of effective (local) exponents df^{L) that describe the 
scaling oi (l) within intervals of, say, 3 successive values 
of L. The change of the effective exponents for increas- 
ing system sizes further show how the scaling behavior is 
affected by the finite size of the simulated systems. From 
the sequence of effective exponents one can extrapolate 
the asymptotic fractal dimension by means of a straight 
line fit to the plot of df^{L) against the inverse system 
size 1/L. Figure [5] shows the effective exponents obtained 
for 3 and 4 successive values of L at different values of 
the disorder parameter p. Therein, the asymptotic frac- 
tal dimensions d{, as listed in table HI where estimated 
from the effective exponents resulting from intervals of 
4 successive system sizes. The asymptotic values for d^, 
9i and 62, listed in tables H^jlll where estimated using a 
similar procedure. Our results for the fractal dimension 
and the stiffness exponent at p = 0.60 and 0.62 clearly 
support the estimates for the pure SG at p = 0. They are 
in agreement with the SLE scaling relation and hence we 
could verify that the SLE scaling relation holds up to 
values of the disorder parameter close to pc- 

At the critical point we find that the estimates of df 
and 62 are not in agreement with the SLE scaling relation. 
However, MEDWs at pc are self-similar with the scaling 




1/L 

FIG. 5: Extrapolation of the asymptotic fractal dimension 
for Model I. Analysis of the sequence of effective exponents 
df^ {L) that describe the scaling of the average MEDW length 
within intervals of 3 (open symbols) and 4 (filled symbols) 
successive values of L, according to {I) ~L'*^ The asymptotic 
value of the fractal dimension di is extrapolated from the plot 
of df * against 1 / L as the intersection of a straight line fit to 
the data with the ordinate. 

dimension df = 1.222(1) and a roughness compatible with 
unity. Here, the numerical value of df as estimated from 
the effective exponents compares nicely to the value df ~ 
1.222(4) found from the previous FSS. 

Wc further find that df and 62 in the ferromagnetic 
phase above the critical point are not consistent with 
the SLE scaling relation. There, the overall length of 
the DW increases linear with the system size, i.e. the 
fractal dimension extrapolates towards df = 1, whereas 
for the roughness exponent dj. < 1 is found. This in- 
dicates that, albeit the MEDW is allowed to bend and 
turn back and forth on the lattice, the resulting over- 
hangs are not significant for their scaling behavior. Fur- 



p 


df 


dr 


6*1 


62 


0.00 


1.274(2) 


1.008(11) 


-0.287(4) 


-0.287(4) 


0.60 


1.275(1) 


1.003(3) 


-0.28(1) 


-0.28(2) 


0.64 


1.275(2) 


1.012(4) 


-0.28(1) 


-0.28(4) 


0.66 


1.222(1) 


1.002(2) 


0.17(2) 


0.16(1) 


0.68 


1.05(2) 


0.74(3) 


0.97(4) 


0.35(3) 


0.72 


1.022(1) 


0.698(6) 


1.052(3) 


0.27(2) 



TABLE L From left to right: disorder parameter, fractal di- 
mension, roughness exponent and exponents that describe the 
scaling of the mean and width of the MEDW energy distri- 
bution. The figures for p = are taken from p^ . 



r 


df 


dr 


61 


02 


00 


1.274(2) 


1.008(11) 


-0.287(4) 


-0.287(4) 


1.111 


1.275(7) 


0.994(4) 


-0.294(6) 


-0.295(5) 


1.010 


1.286(3) 


1.024(2) 


-0.311(2) 


-0.35(1) 


0.970 


1.222(6) 


0.999(3) 


0.15(1) 


0.15(1) 


0.935 


1.085(4) 


0.782(3) 


0.96(2) 


0.31(2) 


0.833 


1.015(1) 


0.651(3) 


1.028(1) 


0.31(2) 



TABLE II: From left to right: disorder parameter r = aj/pj, 
fractal dimension, roughness exponent and exponents that 
describe the scaling of the mean and width of the MEDW 
energy distribution. The figures for r = 00, i.e. pj = ,are 
taken from [l5| . 
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ther, the cost needed to introduce the DW grows almost 
Unearly with the system size, while the rms-fluctuation 
is characterized by an exponent significantly smaller than 
that. Hence, MEDWs in the ferromagnetic phase display 
a self-afEne scaling, governed by exponents that are in 
reasonable agreement with those that describe the scal- 
ing of the transverse deviation {^L'^/^) and the rms ex- 
citation energy {'-^ L^/'^) of pinned DWs in an ordinary 
Ising FM with randomly placed impurities (33| . Further, 
the scaling behavior found here agrees with that observed 
for directed and undirected optimal paths on 2d lattices 
subject to weak disorder (s^l or analogously the scaling 
of directed polymers in random media (35[ . 

For the Model II disorder, our findings are qualitatively 
the same, hence we only state the numerical results with- 
out showing figures. The numerical values of df and 02 
within the SG phase (r > Tc) are in agreement with the 
SLE scaling relation proposed for the pure SG. In partic- 
ular, at r^l.Ol we find df = 1.286(3) and 6*2 = -0.35(1). 
Here, the data for (/) gives a nice straight line on a double 
logarithmic scale, where we find df = 1.284(2) from a fit 
to the pure power law data excluding L < 100. The sit- 
uation for the data corresponding to a{SE) is somewhat 
different, i.e. the data still exhibits a curvature within 
the range of accessible system sizes on a double loga- 
rithmic scale. This does not allow to fit all the data at 
once, assuming a power law fit-function. Consequently, 
the most reliable estimate of the asymptotic value of 6*2 
can be obtained by an analysis of the local exponents 
as described above. The reason for this difficulty might 
stem from the fact that the value r = 1.01 of the disorder 
parameter is located in the transition region close to the 
critical point. Albeit these values differ slightly from the 
values 62 ~ —0.28 and di ~ 1.274 that one would expect to 
find in the SG phase, they are in agreement with the SLE 
scaling relation. The numerical values for the exponents 
right at the critical point are again not in agreement with 
the proposed scaling relation. However, the asymptotic 
fractal dimension extrapolated from the effective expo- 
nents reads df = 1.222(6) and is in agreement with the 
corresponding value at the critical point for the Model I 
disorder. Further, if we analyze the scaling of the average 
MEDW length restricted to system sizes L < 64 we find a 
value of d{ = 1.246(4), consistent with the value encoun- 
tered in the previous FSS analysis that was denoted as a 
finite-size effect. Within the ferromagnetic phase (r < Tc), 
MEDWs again display a self-affine scaling behavior, fur- 
ther characterized by exponents that extrapolate towards 
those that describe the scaling of the transverse devia- 
tion and the rms excitation energy of pinned DWs in an 
ordinary Ising FM with randomly placed impurities. Re- 
ferring to the droplet model and under the assumption 
that all lengths exhibit the same asymptotic scaling be- 
havior, one can further relate the exponents df and 62 
by means of the equation df = d/2 — O2 [1^. Note that 
within the SG phase our data is in reasonable agreement 
with this scaling relation. In this regard, the agreement 
for model I is somewhat better than for model II. Fur- 




FIG. 6: Samples of GS spin configurations for systems of 
side-lengtli L = 64 and r = 1.1, 0.969, 0.9 (from left to right) 
for Model II. 




0.9 0.95 1 1.05 1.1 1.15 



r 

FIG. 7: FSS analysis of ferromagnetic domains at T = 
for Model II. The main plot shows the average size-ratio 
(Ah/Ml) of the second-largest and largest ferromagnetic clus- 
ters and the inset illustrates the data collapse under the re- 
spective scaling assumption, obtained for Vc = 0.969(2) and 
!/F = 1.49(4). 



ther, combining the above equation with the SLE scaling 
relation [sj, we would expect 02 ~ —0.2753 or similarly 
df« 1.2753. 



C. Finite-size scaling analysis of ferromagnetic spin 
domains at the r = spin glass to ferromagnet 
transition 

As we decrease the value of the disorder parameter in 
Model II from r = 00 (SG-phasc) to r < Tc (FM-phase), 
we can identify ferromagnetic clusters of spins, i.e. groups 
of nearest-neighbor spins with similar orientation, with 
increasing size (see Fig. [5]). Here, as an alternative way 
to characterize the SG to FM transition at T = 0, we per- 
form a FSS analysis of the largest and second largest fer- 
romagnetic clusters found for the GS spin configuration 
for each realization of the disorder. Such an analysis has 
been performed previously for standard percolation [ssj . 
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FIG. 8: FSS analysis of ferromagnetic domains at T = 
for Model II. (a) normalized size (Mi) of the largest fer- 
romagnetic cluster, and (b) finite-size susceptibility xl ~ 
N[{Mi) — (A/i)^] associated with the the size of the largest 
cluster. The main plots show the unsealed data near the crit- 
ical point, while the insets illustrate the data collapse under 
the respective scaling assumptions. 



As above, we simulated systems of size L = 24, 32, 48, 64 
at different values of the disorder parameter r. We kept 
the width of the disorder distribution at the fixed value 
aj = I and we vary only its mean fij. Observables are 
averaged over 2x10^ samples for each system size. Subse- 
quently, the relative size of a cluster specifies the number 
of spins that comprise the cluster divided by the num- 
ber of spins on the lattice. Within our analysis we found 
that the average ratio {M2/M1) of the relative sizes of 
the second-largest and the largest ferromagnetic clusters 



scales as 

{M2/M,)^Mr-r,)L'/''-], (6) 

therein Tc is the location of the critical point and 
signifies the correlation length exponent. From a data 
collapse, restricted to the interval [—1.5,4-1.5] on the 
rescaled abscissa, we obtain the numerical values Tc = 
0.969(2) and i^f = 1.49(4) with a quality S = 0.85, see 
figure [T] Both values agree within errorbars with those 
obtained from the Binder parameter analysis. If we allow 
for a nonzero scaling dimension according to {M2I Mi) 
L^'^feKr — rc)L^^'^'^], we yield Tc and lyp as above and 
further k = 0.004(13) (5 = 0.82, [-2.0, +1.0]). The nu- 
merical value of K is compatible with zero and hence sup- 
ports the scaling assumption © for the size ratio. More- 
over, right at Tc we found the critical value {M2/M1) = 
0.122(1). For completeness we note, that we yield qual- 
itatively similar findings for the ratio {M2) / {Mi) with 
the critical value (M2)/(Mi) = 0.098(1). The difference 
between the two ratios is simply due to the cluster-size 
fluctuations at criticality. 

As an order parameter we measure the relative size 
Ml of the largest ferromagnetic cluster for each of the 
GSs. From the scaling assumption {Mi) L-f^^/"^ f^[{r- 
rc)L^^'^^] and the values of Tc and v-p stated above we 
obtain /3f =0.039(4) (5* = 0.54, [-0.5,4-0.5]), see figure 
^i&). A similar scaling assumption for the second largest 
cluster yields /?f,2 = 0.05(3) (5 = 0.26, [-0.5,4-0.25], not 
shown). Albeit the numerical value of /3f,2 is less precise 
and somewhat larger compared to /3f, both exponents 
are compatible with Eq. 

The finite-size susceptibility xl = N[{Ml) - (A/i)^] 
describing the fluctuations of the size of the largest 
ferromagnetic cluster, obeys the scaling form xl ^ 
L^^/'^^ f^[{r — rc)L^^'^^] with another critical exponent 
7f, see flgure[8Ub). Together with the values of Tc and v-p 
we estimate 7F = 2.9(1) (5 = 0.85, [-1.5, +1.0]). These 
exponents further are in agreement with the hyperscaling 
relation 7f/^'f 4- 2Pp/v-p = d. 

We performed further simulations for the ±J model 
with a varying fraction 0.0 < p < 0.5 of aniferromagnetic 
bonds. In principle, the GS for this model is highly de- 
generate [H,!!^]. Here, we investigate only one randomly 
obtained GS for each realization of the disorder. From a 
FSS analysis for systems of size L = 32, 48, 64, 96, where 
averages are computed over 3x10^ samples, we found Pc = 
0.1022(3), i^F = 1.47(6), (ip = 0.037(4) and 7f = 2.8(1). 
The numerical values of the critical exponents for the 
±J model agree, within errorbars, with those obtained 
for Model II above. Further, the critical concentration 
of antiferromagnetic bonds is in fair agreement with the 
value Pc = 0.103(1) found from an analysis of the Binder 
parameter within a previous study |3l| . Regarding the 
FSS analysis and compared to [3l| , we used a larger num- 
ber of interpolation points that enclose the critical point 
on the rescaled abscissa (24 data points in the interval 
[—0.5 : 0.5] for each system size). As a result we obtained 
Pc with increased precision, although our system sizes are 
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somewhat smaller. Finally, right at Pc we found the criti- 
cal ratios (M2/Mi) = 0.104(1) and (M2)/(Mi) = 0.083(1). 
The numerical values of these ratios differ slightly from 
those obtained for Model II above. However, in both 
cases we observe (M2/Afi) «0.125 {M2)/{Mi). 

As mentioned above, the scaling of the size ratio ac- 
cording to equati on El was also confirmed for usual ran- 
dom percolation [S^- It stems from the fact that the 
largest and second-largest clusters exhibit the same frac- 
tal dimension at the critical point. For usual percolation 
this was shown earlier [4l|. While we could verify equa- 
tion © for the disorder induced SG to FM transition 
at T = numerically, we found within additional simu- 
lations no such scaling behavior for the thermal phase 
transition in the 2d Ising ferromagnet. 



IV. SUMMARY 

We have investigated MEDWs for two-dimensional 
random-bond Ising spin systems, regarding two differ- 
ent continuous bond distributions. For both models, a 
disorder parameter could be used to distinguish between 
a spin-glass ordered or a ferromagnetic ground state. We 
performed a FSS analysis to locate the critical points in 
both models that separate the spin-glass phase from the 
ferromagnetic phase. We found that within the spin-glass 
phase, the exponents that describe the size scaling of the 
width of the average DW energy and the average DW 
length are approximately constant and consistent with 
the SLE scaling relation previously proposed for the pure 
spin-glass. Right at the critical point and in the ferro- 
magnetic phase of the models the accordant exponents 
are not in agreement with the SLE scaling relation. 

It is intriguing to note that the fractal dimension of the 
DWs at the critical point of both disorder types studied 
here, agrees with the fractal dimension dopt = 1.22(2) 
of optimal paths in the strong disorder limit on 2d lat- 
tices [42| . This is quite interesting since the optimization 
criteria of the two problems are rather distinct: In the 



strong disorder limit, nonnegative edge weights are drawn 
from a very broad distribution. The cost of a path be- 
tween two sites on the lattice is then dominated by the 
largest edge-weight along the path. Consequently, so as 
to find an optimal path, one has to minimize the largest 



weight along the path. In contrast to this, the cost of a 
MEDW is the sum of all edge weights along the respective 
lattice path. There are positive and also negative edge 
weights that can cancel each other, at least partially. A 
common feature of the above two problems is that, in 
striking contrast to usual shortest path problems, there 
is no immediate negative feedback for the inclusion of 
additional path segments. In usual shortest path prob- 
lems, where there are only positive edge weights, like e.g. 
optimal paths subject to weak disorder [3J], the inclu- 
sion of additional path segments leads very likely to an 
increased path weight. Hence, postive-weight minimum- 
weight paths tend to be short, which results in an average 
end-to-end distance ~L. 

Finally, we have characterized the SG to FM transition 
at T = in terms of the largest and second-largest ferro- 
magnetic clusters of spins found for the GS spin config- 
urations. The respective critical exponents support our 
previous results and they appear to be consistent with a 
hyperscaling relation known from scaling theory. 

After submission of our manuscript to Phys. Rev. B, 
we received correspondence from R. Fisch and E. Vicari 
containing valuable suggestions and comments that en- 
abled us to improve the manuscript further. The article 
was accepted without a request for amendments. Hence 
and to our regret, these improvements do not appear in 
the published version [Phys. Rev. B 79, 184402 (2009)] 
of the article. 
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